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The three dimensional (3D) topological insulators are predicted to exhibit a 3D Dirac semimetal 
state in critical regime of topological to trivial phase transition. Here we demonstrate the first 
experimental evidence of 3D Dirac semimetal state in topological insulator Bi 2 Se 3 with bulk carrier 
concentration of 10^^ cm“®, using magneto-transport measurements. At low temperatures, the 
resistivity of our Bi 2 Se 3 crystal exhibits clear Shubnikov-de Haas (SdH) oscillations above 6 T. The 
analysis of these oscillations through Lifshitz-Onsanger and Lifshitz-Kosevich theory reveals a non¬ 
trivial TT Berry phase coming from 3D bands, which is a decisive signature of 3D Dirac semimetal 
state. The large value of Dingle temperature and natural selenium vacancies in our crystal suggest 
that the observed 3D Dirac semimetal state is an outcome of enhanced strain field and weaker 
effective spin-orbit coupling. 


The 3D topological insulators are an exotic class of 
quantum materials having insulating bulk and novel 
spin-momentum locked dissipationless conducting sur¬ 
face. The surface states of these materials consist of odd 
numbers of massless Dirac fermions with linear disper¬ 
sion protected by time reversal symmetry 0 , 0 . In 3D 
topological insulators, the chalcogenide 612803 has drawn 
ample attention because of relatively simple and robust 
surface states consisting of single Dirac cone at T point of 
band structure and a large topological band gap (0.3 eV) 
between the bulk states 0. The existence of topologi¬ 
cal surface states in 612803 has been conhrmed through 
angle resolved photoemission spectroscopy (ARPE8) 0 , 
scanning tunneling microscopy 0, and transport mea¬ 
surements 0. The as grown crystals of 612803 are typ¬ 
ically n-type because of electron doping due to natural 
selenium vacancies 0[3|. Therefore, the transport prop¬ 
erties of 612803 are generally dominated by bulk conduc¬ 
tion, for example, the temperature dependence of resis¬ 
tivity is metal like [M and 8hubnikov-de Haas (SdH) 
oscillations in resistivity have the characteristic of bulk 
Fermi surface ISHi. The bulk carrier concentration (n) 
and Fermi energy in 612803 can be tuned by either con¬ 
trolling the natural Se vacancies 0 0 , gate-voltage 0, 
or chemical doping (bl-lsl. [l2|. For n^lO^^cm SdH os¬ 


cillations are from two dimensional (2D) states indicating 
the surface transport 0 , while for n>10^^cm“^, SdH os¬ 
cillations come from 3D states indicating the dominance 
of bulk conduction 01i3|. For high n i.e. n~10^®cm 
bulk dominated transport with 2D like SdH and quantum 
Hall effect is reported [10. 

The Dirac fermions at the surface of 3D-topological 
insulators are two dimensional i.e. graphene like, and 
very recently, their three dimensional analog have been 
discovered in a new kind of Dirac material known as 3D 
Dirac semimetal. In 3D Dirac semimetals, the bulk con¬ 
duction and valance band touch at the discrete Dirac 
points and disperse linearly along all momentum direc¬ 
tions The 3D Dirac semimetals are predicted 

to exist at the critical point of topological phase transi¬ 
tion from topological to trivial -insulator in a space in¬ 


version symmetry 0. 8uch a topological phase tran¬ 
sition can be realized by closing the bulk band gap at 
the critical point by tuning of effective spin-orbit cou¬ 
pling with chemical doping [isl - l^ or varying the lat¬ 
tice parameters through strain/pressure [^-24|. Exper¬ 
imentally the 3D Dirac semimetal state has been ob¬ 
served by ARPE8 measurements on EiTl(8o.5800.5)2 0, 
(Ein Q^Inn nfi)98 e3 [l9l| . and in T16i88e through transport 
measurements [20. 

In this letter we report the quantum transport evidence 
of 3D Dirac semimetal state in as prepared Ei2Se3 single 
crystal. We show that the SdH oscillations in Bi2Se3 
crystal having carrier concentration of ^1.13x 10^® cm“^, 
originate from 3D bands and has a non-trivial tt Berry 
phase. The presence of tt Berry phase from 3D bands 
suggest that the electrons encircle the 3D Dirac point 
in the momentum space. The large Dingle temperature 
in our crystal indicates the presence of high strain field. 
This transformation from topological insulator to Dirac 
semimetal state is a result of enhanced strain field and 
weaker effective spin-orbit coupling. 

Bi28e3 single crystal isgrown by modified Bridgeman 
method as in Reference [ISj . The as grown crystal is eas¬ 
ily cleaved which yields a flat shiny surface. The crystal 
is characterized by X-ray diffraction (XRD) performed on 
a Bruker D8 Advance X-ray diffractometer using Cu Ka 
radiation, the data is shown in Figured] (a). The XRD 
data exhibit peaks only corresponding to (003) family 
of planes, suggesting the cleaved surface is perpendicu¬ 
lar to C3 axis of rhombohedral structure. The powder 
X-ray diffraction of the crushed crystal is displayed in 
Figured) (b)- The Le Bail Fitting of the powder XRD 
data gives a=4.144 A and c=28.660 A. The energy dis¬ 
persive X-ray analysis (EDX) yields atomic ratio of Bi 
and 8e as 2 and 2.84 respectively. 

A thin rectangular bar shaped Bi2Se3 crystal was taken 
and contacts for resistivity and Hall measurements were 
made on freshly cleaved surface. Resistivity measure¬ 
ments are performed using standard four probe ac tech¬ 
nique while Hall measurements are carried through hve 
probe ac technique on a 9T PPM8 AC Transport system 
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Figure 1: X-ray diffraction pattern of (a) freshly cleaved 
Bi 2 Se 3 single crystal and (b) Bi 2 Se 3 powder obtained 
after crushing the crystal. The symbols correspond to 
experimental data, red line represents the Le Bail fitting 
{R3fh space group) and blue line is the difference 
between fitted and experimental data. 



Figure 2: (a) Temperature dependence of longitudinal 
resistivity (pxx) m zero magnetic field, (b) pxx and (c) 
Hall resistance {pxy) as a function of magnetic field 
(HIIC 3 axis) at 1.8 K. (d) Expanded view of pxx versus 
B plot exhibiting clear Shubinikov-de Haas oscillations. 


(Quantum Design). The angle dependent resistivity mea¬ 
surements are performed on the high resolution rotator 
option of PPMS. 

Figure [2] (a) shows the temperature dependence of re¬ 
sistivity Pxx- The temperature variation of resistivity 
displays metallic character with residual resistivity ratio 
(RRR) of about 2.7 which lies in the range of ~1.6-3.0 
as observed in earlier reports [3, [ll|. Figure [11(b) and 
(c) display the field dependence (field parallel to C 3 axis) 
of longitudinal resistivity (pxx) and Hall resistivity (pxy) 
at 1.8 K respectively. The longitudinal resistivity {pxx) 
increases on increasing the magnetic field, and above 6 T, 
exhibits periodic oscillations with field. An expanded 
view of these oscillations is shown in figured] (d). The 



Figure 3: Oscillatory component of longitudinal 
resistivity (Apxx) as a function of magnetic filed at 
different temperatures. Apxx is obtained after 
subtracting a smooth background from pxx- 


resistivity oscillations in metals occur due to formation 
of Landau levels at high magnetic fields and are known as 
Shubnikov-de Haas (SdH) oscillations. The result of field 
variation of Hall resistivity shows that our sample is n- 
type and a linear fitting of the data above 4 T gives bulk 
carrier density n|^“=1.13x 10^® cm“^. The Hall mobil¬ 
ity (ph) estimated from ratio of Hall coefficient and lon¬ 
gitudinal resistivity is ^986 cm^Vs“^. The n-type doping 
in Bi 2 Se 3 occurs due to natural Se vacancies. 

Figure |3| displays the SdH oscillations in Ap {Ap is 
obtained after subtracting a polynomial background fit 
from Pxx) above 6 T in the temperature range of 1.8 
to 12 K. The amplitude of SdH oscillations decreases 
on increasing the temperature and the oscillations are 
not resolvable above 12 K. The oscillations have a sin¬ 
gle frequency Bp=219{l) T which is large in comparison 
to the previously reported values for Bi 2 Se 3 single crys¬ 
tals [^, [ 1 ^, El and nano-flakes (26l - [2^ . The Landau level 
index n of the SdH oscillations is related to extremal cross 
section area of Fermi surface (T^kp) perpendicular to the 
applied field through Lifshitz-Onsanger equation 

2Tr{n + j)=Trkj,-^, ( 1 ) 

where h is the reduced Plank’s constant, e is the elec¬ 
tronic charge, and B is the applied magnetic field. The 
coefficient 7, which is equal but opposite in sign to the 
intercept on Y-axis in n versus 1/B plot, is related to the 
phase of SdH oscillations. In absence of iiitraband and 
interband magnetic breakdown, for conventional metals 
7=1/2 [^. 7 gets an additional non-trivial contribution 
from Berry phase when electron orbit in k space sur¬ 
rounds a contact line of its band with some other band, 
and with Dirac like dispersion, these band contact lines 
form the Dirac points |30|. In general, the phase of SdH 
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Figure 4: (a) Landau index n versus 1 /B plot for Ap at 
1.8 K. The straight line is the least square fitting of 
Lifshitz-Onsanger equation giving 7=0.15(5). Inset 
shows the magnified view around the intercept, (b) 
Direct fitting of Lifshitz-Kosevich equation to Ap at 
1.8 K. 


oscillation 7 = 1/2 — $B/ 27 r — 5 , where is the Berry 
phase (<i>B=0 for conventional metals and tt for Dirac sys¬ 
tems) and 2tt5 is the phase shift correction determined by 
the dimensionality of Fermi surface ^ l32|. 5=0 for 2D 
Fermi surface and 5=±l/8 for 3D Fermi surface with pos¬ 
itive sign for holes and negative for electrons. Figure|4](a) 
exhibits the Landau level fan diagram, i.e. Landau level 
index n versus 1/B plot, for SdH oscillations at 1.8 K. 
The minima of Ap versus B plot are assigned integer val¬ 
ues of n while half integer values n+ 1/2 are assigned to 
maxima (s^, [s^ . The data in figure 0 (a) fits well with 
the equation [T] giving 7=0.15(5) and /cf=0.0815(1) A“^. 
7 should be 1/2 for conventional metals, 0 for 2D and 
^1/8 (0.125) for 3D Dirac electrons with tt Berry phase. 
Therefore, 7=0.15(5) gives strong evidence for the pres¬ 
ence of non-trivial tt Berry phase arising from 3D Dirac 
electrons. The 7 value is further reassessed by fitting 
the Ap oscillations directly with Lifshitz-Kosevich equa¬ 
tion [3i,[3i| 


Ap{B) = Apo{T)exp{—TT / pB)cos2Tr{Bp / B + 7). (2) 

The htting of equation to SdH oscillations at 1.8 K 
is shown in Figure S] (b). We obtain htting parameters 


Figure 5: (a) SdH oscillations in Ap at various tilt 
angles (0) at 2 K. (b) The frequency of SdH oscillations 
Bi? as a function of tilt angle (Black squares). The red 
line with circles is the expected angle dependence of B^? 
for 2D Fermi surface (B/^’ oc l/cos0). The blue solid line 
is the least square htting of prolate spheroid Fermi 
surface to B^ versus 9. 


7 = 0 . 2 ( 1 ), Bf=220{1), and p=8{l)xl0^ cm^/Vs. The 
7 value from direct htting reconhrms the result of Lan¬ 
dau level fan diagram. SdH oscillations give carrier con¬ 
centration of 5.3x10^^ cm“^ for 2D Fermi surface and 
1.8x10^® cm“^ for 3D Fermi surface. The mobility and 
carrier concentration (for 3D Fermi surface) values ob¬ 
tained from SdH oscillations are close to their corre¬ 
sponding values from Hall measurements which further 
support the origin of SdH oscillations from 3D bands. 
The relatively large value oi Bp (and kp) in our sys¬ 
tem indicates that the SdH oscillations are coming from 
a large Fermi surface. 

The angle dependence of the SdH oscillations is used 
to determine the shape and size of the Fermi surface. 
Figure [S] (a) shows the SdH oscillations at different tilt 
angles ( 0 ) while hgure [S] (b) exhibits the angle depen¬ 
dence of SdH oscillation frequency at 2 K. Here 6 is 
the angle between C 3 axis and magnetic held; 6 is zero 
when magnetic held is along the C 3 axis and 90° when 
magnetic held is perpendicular to the C 3 axis. For 2D 
Fermi surface, the SdH oscillation frequency oc l/cos0. 
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Figure 6: (a) The temperature dependence of SdH 
oscillation amplitude Ap for n=25+l/2. The solid line 
is a fit of equation [ 3 ] which gives m*=0.15me. (b) 
Dingle plot of \nD versus 1/ B at 1.8 K. 


The data shown in hgure [5] (b) strongly deviates from 
l/cos0 dependence which confirms that the SdH oscil¬ 
lations in our system are not from 2D surface states. 
The angle dependence of SdH oscillation frequency fits 
reasonably well with prolate spheroid Fermi surface 

Bp={h^ky/2e)^J {kl/ky) sin^0 + cos^O with major axis 

A:^=0.112(2) and minor axis kx=ky=0.Q'&1{2) A“^. 
These values of kx and ky are almost same as that re¬ 
ported for topological surface states in ARPES measure¬ 
ments (fcj;=0.09 A“^ and ^2^=0.10 A“^) [3|. The angle 
dependence of SdH oscillations also confirms the origin 
of SdH oscillations in our system is from 3D bands. 

Figure [ 3 ] (a) displays the temperature dependence 
of the amplitude of SdH oscillation for n=25-|-l/2 
(B«8.85 T). The amplitude of SdH oscillation reduces on 
increasing the temperature and totally diminishes above 
12 K. According to Lifshitz-Kosevich theory, the temper¬ 
ature damping of the amplitude of SdH oscillation can 
be given as [sj, 1 ^ : 


Apo(T) = Apo(O) 


27r^fcBT /hujc 
sinh(27r2fcsT/ha;c) ’ 


( 3 ) 


where Apo(O) is the oscillation amplitude at T=0, kp 
is the Boltzman constant, and uJc is the cyclotron fre¬ 
quency. uJc=eB/m* where m* is the electron effective 
mass. The fitting of equation|3]to the data of Figure|n](a) 
yields TO*=0.15me, where me (=9.lx 10“^^ kg) is the free 
electron mass. The value of electron effective mass is 
in agrement with the previous reports (1, [l^, [3|- The 
energy difference between two consecutive Landau levels 
AE=huJc 6.8 meV for 8.85 T. Using the value of elec¬ 
tron effective mass (m*) and Fermi wave vector (kp), the 
Fermi velocity vp=hkp/m* is estimated as 6.3x10® m/s. 
The value of Fermi velocity is similar to that reported for 
surface states in ARPES measurement (« 5 x 10®m/s) [3|- 
The Fermi energy Ep determined from relativistic for¬ 
mula Ep=m*Vp ^ 336 meV above the Dirac point. 

The Dingle temperature (To) which accounts for the 


broadening of Landau levels due to electron scattering 
is given as Tp)= h/2'KkBTs, where Ts is the quantum life 
time of carrier due to scattering. Figure[n](b) displays the 
semilog plot of D=ApBsmh{27r'^kBT/fiuJc) versus 1/B at 
1.8 K and the Dingle temperature is estimated from the 
slope of linear fit which is ~40 K. From Dingle tempera¬ 
ture we get the carrier life time rs?a3x 10“^'’’ s. The Din¬ 
gle temperature of our sample is large in comparison to 
the previously reported values for Bi 2 Se 3 e.g. T£)=3.5 K 
for na!2.3xl0i^ cm"® Td= 9.5 K for n^SxlO^® 
cm“®[^, Tp)=4.0 K for nRi2.3xlO^® cm“® and 
T£)=25 K for nwSxlO^® cm“®[^. The large value of 
T B (for similar carrier concentration) in our sample is 
an indication of the presence of relatively high strain 
fields (^. The strain field plays a vital role in determin¬ 
ing the topological properties of Bi 2 Se 3 , and it has been 
observed that the strain field may drive the Bi 2 Se 3 sys¬ 
tem through topological to trivial -insulator phase tran¬ 
sition [i^ . 


The observation of tt Berry phase in 3D bands confirms 
the presence of Dirac point with 3D electron dispersion. 
This suggests that our Bi 2 Se 3 system is in a 3D Dirac 
semimetal state. The 3D Dirac semimetal state in Bi 2 Se 3 
can be realized on closure of the bulk topological band 
gap in the critical regime of topological phase transition. 
It has been shown that the bulk topological band gap 
in Bi 2 Se 3 can be reduced on lowering the effective spin- 
orbit coupling as well as on enhancing the strain 

field in the crystal 1J,12^. In our system, atomic ratio of 
Bi:Se is 2:2.84 due to natural Se vacancies, suggesting the 
possibility of a slightly lower effective spin-orbit coupling 
than that of perfect stoichiometric crystal. Furthermore, 
the observation of large Dingle temperature in compari¬ 
son to pervious reports asserts the presence of enhanced 
strain field in the system. These results imply that the 
3D Dirac semimetal state in our system is an outcome of 
the closing of bulk topological band gap due to weaker 
spin-orbit coupling and enhanced strain field. 


In summary, magneto-transport measurements on our 
Bi 2 Se 3 single crystal give strong evidence of non-trivial tt 
Berry phase originating from 3D bands. The observation 
of non-trivial tt Berry phase from 3D bands suggests that 
our Bi 2 Se 3 system is in a 3D Dirac semimetal state. The 
analysis of Dingle temperature and Bi:Se atomic ratio in¬ 
dicate that the transformation from topological insulator 
to Dirac semimetal state is because of enhanced strain 
field and reduced effective spin-orbit coupling. The rev¬ 
elation of Dirac semimetal state would provide a new in¬ 
sight in understanding the electronic properties of topo¬ 
logical insulators and the tunability of topological phase 
without external doping, which preserves the high elec¬ 
tron mobility, would offer interesting opportunities for 
future device applications. 

















5 


ACKNOWLEDGEMENTS 

We thank A. Tayal and M. Gupta for XRD measure¬ 
ments and D. M. Phase for EDX measurements. V. 
Ganesan and A. Banerjee are acknowledged for support. 


* Electronic address: deveniit@gmail.com 

[1] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

( 2010 ). 

[2] X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 

( 2011 ). 

[3] H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. 
Zhang, Nature Physics 5, 438 (2009). 

[4] Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. 
Bansil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. 
Hasan, Nature Physics 5, 398 (2009). 

[5] P. Cheng, C. Song, T. Zhang, Y. Zhang, Y. Wang, J.-F. 
Jia, J. Wang, Y. Wang, B.-F. Zhu, X. Chen, X. Ma, K. 
He, L. Wang, X. Dai, Z. Fang, X. Xie, X.-L. Qi, C.-X. 
Liu, S.-C. Zhang, and Q.-K. Xue, Phys. Rev. Lett. 105 , 
076801 (2010). 

[6] J. G. Analytis, R. D. McDonald, S. C. Riggs, J.-H. Chu, 
G. S. Boebinger, and Ian R. Fisher, Nature Physics 6, 
960 (2010). 

[7] Y. S. Hor, A. Richardella, P. Roushan, Y. Xia, J. G. 
Gheckelsky, A. Yazdani, M. Z. Hasan, N. P. Ong, and R. 
J. Cava, Phys. Rev. B 79 , 195208 (2009). 

[8] J. G. Analytis, J.-H. Chu, Y. Chen, F. Corredor, R. D. 
McDonald, Z. X. Shen, and 1. R. Fisher, Phys. Rev. B 
81, 205407 (2010). 

[9] N. P. Butch, K. Kirshenbaum, P. Syers, A. B. Sushkov, 
G. S. Jenkins, H. D. Drew, and J. Paglione, Phys. Rev. 
B 81 , 241301(R)(2010). 

[10] K. Eto, Z. Ren, A. A. Taskin, K. Segawa, and Y. Ando, 
Phys. Rev. B 81, 195309 (2010). 

[11] J. G. Gheckelsky, Y. S. Hor, R. J. Cava, and N. P. Ong, 
Phys. Rev. Lett. 106 , 196801 (2011). 

[12] Z. Ren, A. A. Taskin, S. Sasaki, K. Segawa, and Y. Ando, 
Phys. Rev. B 84, 075316 (2011). 

[13] H. Cao, J. Tian, I. Miotkowski, T. Shen, J. Hu, S. Qiao, 
and Y. P. Chen, Phys. Rev. Lett. 108, 216803 (2012). 

[14] S. M. Young, S. Zaheer, J. C. Y. Teo, C. L. Kane, E. J. 
Mele, and A. M. Rappe, Phys. Rev. Lett. 108 , 140405 
( 2012 ). 

[15] B.-J. Yang and N. Nagaosa, Nature Communications 5, 
4898 (2014). 


[16] L. P. He, X. C. Hong, J. K. Dong, J. Pan, Z. Zhang, 
J. Zhang, and S. Y. Li, Phys. Rev. Lett. 113, 246402 
(2014). 

[17] S. Murakami, New J. Phys. 9, 356 (2007). 

[18] S.-Y. Xu, Y. Xia, L. A. Wray, S. Jia, F. Meier, J. H. Dil, 

J. Osterwalder, B. Slomski, A. Bansil, H. Lin, R. J. Cava, 
M. Z. Hasan, Science 332, 560 (2011). 

[19] M. Brahlek, N. Bansal, N. Koirala, S.-Y. Xu, M. Neu- 
pane, C. Liu, M. Z. Hasan, and S. Oh, Phys. Rev. Lett. 
109, 186403 (2012). 

[20] J. Liu and D. Vanderbilt, Phys. Rev. B 88, 224202 
(2013). 

[21] T. Sato, K. Segawa, K. Kosaka, S. Souma, K. Nakayama, 

K. Eto, T. Minami, Y. Ando, and T. Takahashi, Nature 
Physics 7 , 840 (2011). 

[22] Y. Liu, Y. Y. Li, S. Rajput, D. Gilks, L. Lari, P. L. 
Galindo, M.Weinert, V. K. Lazarov, and L. Li, Nature 
Physics 10, 294 (2014). 

[23] W. Liu, X. Peng, G. Tang, L. Sun, K. Zhang, and J. 
Zhong, Phys. Rev. B 84, 245105 (2011). 

[24] S. M. Young, S. Chowdhury, E. J. Walter, E. J. Mele, 

C. L. Kane, and A. M. Rappe, Phys. Rev. B 84, 085106 

( 2011 ). 

[25] M. Novak, S. Sasaki, K. Segawa, and Y. Andov, Phys. 
Rev. B 91, 041203(R) (2015). 

[26] Y. Yan, L.-X. Wang, X. Ke, G. V. Tendeloo, X.-S. Wu, 

D. -P. Yu, Z.-M. Liao, Scientific Reports 4, 3817 (2014). 

[27] Y. Yan, Z.-M. Liao, Y-B. Zhou, H.-C. Wu, Y.-Q. Bie, 
J.-J. Chen, J. Meng, X.-S. Wu, and D.-P. Yu, Scientific 
Reports 3, 1264 (2013). 

[28] D. Kong, J. J. Cha, K. Lai, H. Peng, J. G. Analytis, S. 
Meister, Y. Chen, H.-J. Zhang, 1. R. Fisher, Z.-X. Shen, 
and Yi Cui, ACS Nano 5, 4698 (2011). 

[29] D. Shoenberg, Magnetic Oscillations in Metals (Cam¬ 
bridge University Press, Cambridge, England, 1984). 

[30] G. P. Mikitik and Yu.V. Sharlai, Phys. Rev. Lett. 82, 
2147 (1999). 

[31] H. Murakawa, M. S. Bahramy, M. Tokunaga, Y. Kohama, 
C. Bell, Y. Kaneko, N. Nagaosa, H. Y. Hwang, Y. Tokura, 
Science 342, 1490 (2013). 

[32] Y. Zhao, H. Liu, C. Zhang, H. Wang, J. Wang, Z. Lin, 
Y. Xing, H. Lu, J. Liu, Y. Wang, S. Jia, X. C. Xie, J. 
Wang, arXiv:1412.0330 

[33] D.-X. Qu, Y. S. Hor, J. Xiong, R. J. Cava, N. P. Ong, 
Science 329, 821 (2010). 

[34] S. Barua, K. P. Rajeev and A. K. Gupta, J. Phys.: Con- 
dens. Matter 27 , 015601 (2015). 

[35] H. Tang, D. Liang, R. L. J. Qiu, and X. P. A. Gao, AGS 
Nano 5, 7510 (2011). 



